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ABSTRACT 

Primordial stars are likely to be very massive >30Mq, form in isolation, and will 
likely leave black holes as remnants in the centers of their host dark matter halos. 
We expect primordial stars to form in halos in the mass range 10 6 — 10 10 M©. Some 
of these early black holes, formed at redshifts z>10, could be the seed black hole for 
a significant fraction of the supermassive black holes found in galaxies in the local 
universe. If the black hole descendants of the primordial stars exist, their mergers 
with nearby supermassive black holes may be a prime candidate for long wavelength 
gravitational wave detectors. We simulate formation and evolution of dark matter halos 
in ACDM universe. We seed high-redshift dark matter halos with early black holes, 
and explore the merger history of the host halos and the implications of black hole's 
kick velocities arising from their coalescence. The central concentration of low mass 
early black holes in present day galaxies is reduced if they experience even moderate 
kicks of tens of kms -1 . Even such modest kicks allow the black holes to leave their 
parent halo, which consequently leads to dynamical friction being less effective on the 
low mass black holes that were ejected, compared to those still embedded in their 
parent halos. Therefore, merger rates with central supermassive black holes in the 
largest halos may be reduced by more than an order of magnitude. Using analytical 
and illustrative cosmological N-body simulations, we quantify the role of kicks on the 
merger rates of black holes formed from massive metal free stars with supermassive 
black holes in present day galaxies. 
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1 MOTIVATION 

It is firmly established that most galaxies have super-massive black holes (SMBH) at their centers. Their observed masses 
are in the range 10 6 M <, M £ 10 9 M Q and it appears that there are correlations (Kormendy & Gebhardt 2001, Gerhard 
2001) between their masses and the bulk properties of the galactic bulges hosting them (Laor 2001, Merritt & Ferrarese 2001, 
Tremaine et al. 2002, Gebhardt et al. 2001, Graham et al. 2001). These correlations point to a link between the formation 
of SMBHs and the evolution of their hosts. It also appears that SMBHs are linked to the properties of the host dark matter 
halos. If the SMBH precursors were present from very early on, then their mergers, together with growth by accretion, could 
account for the abundance of the SMBHs today (Schneider et al. 2002). 

Ab initio numerical simulations of the formation of the first luminous objects within the standard structure formation 
framework, finds that metal free stars form in isolation, and may have masses 30 Mq ^ m ^ 300 Mq (Abel et al. 2000, Abel 
et al. 2002, Bromm et al. 2002). In current models of structure formation, dark matter initially dominates and pregalactic 
objects form from small initial density perturbations. As they assemble via hierarchical merging, the metal-free primordial 
gas cools slowly through the rotational lines of hydrogen molecules. As the gas cools, it sinks to the center of the dark matter 
potential well. However, as the cooling rates are small, the clumps fail to fragment as they collapse. This leads to a top-heavy 
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initial stellar mass function (IMF) and to the production of very massive stars, unlike the modern stellar IMF which declines 
rapidly with increasing mass. 

Metal-free, high mass Population III stars lose only a small fraction of their mass in their lifetime of ~3 Myr. This 
suggests that a large population of primordial massive black holes (MBH, with mass <J 10 3 M Q ) could be an end product 
of such pregalactic star formation (Heger et al. 2003). Since they form in rare high-cr density peaks (Couchman & Rees 
1986, Madau &: Rees 2001), relic MBHs are predicted to cluster in the cores of more massive halos (Abel et al. 2002) which 
in turn are formed by subsequent mergers of the original parent dark matter halos. It has been suggested (Volonteri et al. 
2003, Islam et al. 2003) that mergers of massive halos and clustering of these MBHs should start at redshifts as large as 
z~20. Subsequent growth of MBHs would proceed through accretion of baryonic gas and coalescence, possibly leading to the 
formation of intermediate mass black holes (IMBH, with mass in range 10 3 Mq ^ m ^ 10 6 Mq). In the most optimistic 
scenario (Haiman & Loeb 2001), 10 9 Mq black holes could form as early as z~10. 

In this paper we will explore the evolution of the seed black hole in dark matter halos. Mergers of dark matter halos hosting 
at least one black hole (Bromm & Loeb 2003) at high redshift form black hole binaries, and their subsequent coalescence 
under gravitational radiation may give them a significant kick velocity (Favata et al. 2004, Merritt et al. 2004). These recoils 
are caused by anisotropic emission of gravitational waves which carry away linear momentum. The recoil velocity is highly 
controversial requiring full general relativistic calculations. However, recent estimates place most recoils in the range 10-100 
kms -1 , although kicks of a few hundred kms -1 are not unexpected, and the largest recoils should not be above 500 kms -1 
(Favata et al. 2004, Blanchet et al. 2005, see also Baker et al 2006). The amplitude of the kick determines if the black hole 
will be ejected from its host dark matter halo. Even the most massive dark matter halo at redshift z>ll can not retain a 
black hole that receives > 150km s _1 kick, while at redshift of z>8, kicks of 300km s -1 are required to cause ejection (Merritt 
ct al. 2004, Fig 3.). This is valid for black holes with masses 10 3 M <, m <, 10 8 Mq. Kicks of 40km s" 1 should be sufficient 
to: eject IMBHs from globular clusters; displace IMBHs from the centers of dwarf galaxies; and perhaps provide a sufficient 
population of IMBHs for merging at the centers of spiral galaxies (Merritt et al. 2004). 

We use N-body TREEcode simulations (GADGET, Springel et al. 2001) to investigate the observational implications of 
black hole kicks through the asymmetric emission of gravitational radiation during the merger. In this academic exercise, we 
assume that Population III stars collapse to form massive black holes at high redshifts. These black holes grow through gas 
accretion and mergers at the centers of dark matter halos. We use collisionless particles at the centers of dark matter halos 
as tracers of IMBHs. If kicks are important, the growth of IMBHs into supermassive black hole will be suppressed. We test 
two extreme scenarios for this behaviour and whether there are predictable consequences. In the first scenario, kicks do not 
occur in any dark matter halo at any redshift. In the second, kicks occur in every dark matter halo at an initial redshift 
and we introduce a prescription for assigning kicks to our tracer particles. As a result, these two scenarios set limits to the 
IMBHs' merger rate. The first scenario provides a maximum merger rate, and the second scenario provides a lower bound 
on the merger rate (the absolute lower bound is that no mergers occur, for example if black hole formation is suppressed 
and primordial stars leave no compact remnants). A better approach would be to track every individual merger, and assign 
kicks when mergers occur. Our approach is computationally less expensive and provides only constraints. We have chosen 
z=8 as an initial redshift for assigning kicks to our tracers for IMBHs. At redshift z~12 supernovae rates for Pop III stars 
fall to zero (Wise & Abel 2005). The massive black holes resulting from Pop III collapse have already settled in the centers 
of dark matter halos by redshift z~8, the initial redshift for kicks. As mentioned before we are calculating constraints for the 
merger rates and not studying individual mergers. The MBH tracers are tracked from the initial to the final redshift and used 
to quantify the changes in distributions of IMBHs in these two scenarios. Through detailed analysis of their trajectories, we 
study how IMBHs respond to the dynamical evolution of dark matter halos, testing extreme ranges of outcome to check for 
the plausible scope for detection. Our results also apply to black holes ejected by the gravitational slingshot mechanism. The 
slingshot velocities fall in the recoil range (Aarseth 2005). 

Our numerical simulations are followed by an analytical treatment of dynamical friction. As long as the dark matter halos 
have a bound core with density greater than the local galactic density and have "many" particles in them, the numerical 
dynamical friction will mimic "real" dynamical friction with reasonable fidelity (Vine & Sigurdsson 1998). If the halo density 
becomes smaller than the local galactic density, or if the secondary halo structure is down to being resolved by only "few" 
bound particles, then the dynamical friction prescription breaks down and relaxation due to numerical noise and large scale 
potential fluctuations dominate; integration of the dynamical friction term becomes meaningless. This happens basically where 
we stop the simulations, and that is where we analytically estimate a conservative lifetime to coalescence. 

We are exploring the effects of kicks on IMBH merger rates and the formation of SMBHs that are detected in the nuclei 
of luminous galaxies. Knowledge of IMBH merger rates is essential for the future experiments such as LISA (Danzmann 2003) 
and BBO (Big Bang Observer) (Cornish & Crowder 2005). LISA will be able to study much of the last year of inspiral of 
IMBH into SMBHs, as well as the waves from the final collision and coalescence of IMBH binaries. For binary black holes with 
masses in the range 100 M <, M <, 10 4 M Q , LISA can observe the last few years of inspiral, but not the final collision even 
at cosmological distances. Equal-mass black hole binaries enter LISA's frequency band roughly 1000 years before their final 
coalescence, independent of their masses, for the range 100 Mq <5 M <5 10 6 M (Thorne 1995, Cornish & Levin 2002, Vecchio 
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et al. 2004). Third generation gravitational wave observatories, such as BBO mission will have arcsecond-arcminute precision 
in identifying and localizing every merging black hole binary independent of redshift, anywhere in the universe over the years 
to months before their actual merger. In combination with ground-based detectors at higher frequencies, BBO could measure 
the mass, angular momentum, and dynamic spacetime structure of these compact objects with unprecedented precision. 

In §2 we present a straightforward analytical treatment of dynamical friction. In §3 we describe the simulation setup, the 
codes used for analysis, and present results for 3 different high resolution simulations with varying kick parameters for early 
black holes. In §4 we describe the post-merger evolution through dynamical friction based on the arguments in §2, before we 
conclude in §5. 



2 ANALYTICAL EXPECTATIONS 

Clustering of MBHs and their consequent mergers creates a population of IMBHs at the centers of dark matter halos. In this 
section we outline the physical principles acting on a seed black hole over one Hubble time. The early stages of the merger are 
driven by the hierarchical cold collapse of sub-halos (gravitationally bound substructures of larger dark matter halos) within 
a primary halo. We focus our analysis on the single dark matter halo that naturally emerges from the numerical simulations 
as the structure largest in mass in our simulation volume; hence, the primary halo. Subsequent dynamical evolution of the 
IMBH population occurs through dynamical friction, and secular orbital evolution in the presence of any residual triaxiality 
after virialization (Madau & Rees 2001, Holley-Bockelmann et al. 2001, Zhao et al 2002). 

We assume that IMBHs embedded in their parent sub-halos experience dynamical friction acting collectively upon the 
entire compact sub-halo containing the IMBH; at least until tidal stripping significantly reduces the effective sub-halo mass 
to that of the central mass only (Weinberg 1989, Vine & Sigurdsson 1998). Tidal dissolution becomes effective only where the 
density of the primary halo is comparable to that of the sub-halo. With characteristic length scales of ~ 100 pc and masses 
of ~ 10 6 M Q , sub-halo densities are of order 10 0±1 M / pc 3 , and the primary halo density reaches such densities only in the 
inner few kpc, at radii r <, 0.01 r„i r (Fig. 6). 

The primary halo is well-approximated by a singular isothermal sphere for radii of order 0.1 r„j r (Madau & Rees 2001), 
the dispersion is near constant and isotropic over a range in radii, and the density profile is close to isothermal, with a steeper 
fall off at larger radii and flattening at smaller radii. 

Using the mean density interior to r, the distance of the sub-halo of mass M from the center of the primary halo, with 
associated circular velocity v c , the timescale for dynamical friction to bring the sub-halo and associated IMBH to the primary 
halo center (Binney & Tremaine 1987, Eq. 7.26) is: 

(1) 

1 InA GM K ' 

In A is the Coulomb logarithm (Binney & Tremaine 1987, table 7.1). 

This is a conservative estimate for the time scale for early stages of dynamical friction; the orbits of the sub-halos are 
not circular and can therefore sink faster than we estimate, and coupling of the dynamical friction to internal degrees of the 
sub-halo also accelerates the evolution, possibly by as much as a factor of two (Weinberg 1989, Vine & Sigurdsson 1998). 
At small radii, mass loss from tidal disruption becomes significant. Nevertheless, we expect the IMBHs to continue sinking 
towards the center of the primary halo, once they get inside r < 0.01 i v i r . Typical mass ratios of early black holes and their 
host halos are > 10~ 4 . 

Consequently for kicks that expel the IMBH from its host-halo at large redshifts, the dynamical friction time scale is at 
least 10 4 times larger since the IMBH is not embedded in a bound sub-halo. From equation (1), we can calculate the radius, 
r s ink, that allows the sub-halo to sink to the center in less than a Hubble time for a given velocity dispersion and sub-halo 
mass. If every sub-halo carries one IMBH at its center, then the merger rate for IMBHs will be a function of the number 
of sub-halos inside r 3 j„fc. If kicks at higher redshifts supply an IMBH with enough velocity to escape from its sub-halo, the 
expected number of IMBHs at r < r sirL k will be smaller, leading to a significant decrease in the IMBH merger rate. We also 
note that, stopping the growth of IMBH via gas accretion (by ejection into lower density regions) changes the X-ray population 
predictions (Agol & Kamionkowski 2002). We illustrate these effects with N-body simulations in the following sections. 



3 ILLUSTRATIVE SIMULATION 
3.1 Simulation Setup 

We simulated a 10 3 Mpc 3 comoving volume in a ACDM universe with Qm~0.3, Qa=0-7 and h=0.7 from redshift z=40 to 
z=l. From initially low resolution simulations we selected a halo with 1.52 x 10 1 M© at redshift z=l. Around that halo we 
selected a sphere of 2 Mpc comoving radius in the initial conditions and reran the simulations from z=40. Using 128 3 on the 
top level and a refinement factor of 4 in the high resolution region, we attain 4.9 xlO 6 high-resolution particles (softening 
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length 2 kpc comoving) for the simulation and 2.0xl0 6 low-resolution particles (softening length 4 kpc comoving) in the rest 
of the box. The mass of each high resolution particle in these simulations is 8.85 x 10 5 M and the mass of each low-resolution 
particle is 5.66 xlO 7 M Q . 

3.2 Analysis Tools 

Every snapshot consists of masses, positions, and velocities of all types of particles at a specific redshift. Plots of their positions 
show clumps of particles that need to be associated with real-world objects. We used the HOP algorithm (Eisenstein & Hut 
1998) to divide the particles into distinct sets such that particles in individual high-density regions are grouped together 
and are separate from those in other regions. HOP assigns density estimates to each particle by using a Gaussian smoothing 
algorithm and then determines which of the particles has the highest density. Each particle is then associated with its highest 
density neighbour and the process continues to higher and higher densities until the highest density particle is reached. All 
particles that hop to the same maximum are placed into the same group. Every particle is assigned to one and only one group. 
We set a density threshold in order to distinguish between the dense halo and its surroundings and we cut out of the group 
all particles with densities less than 5% of the density of the densest particle. We define a halo as a group of more than 50 
particles. The HOP output is a list of halos with IDs, coordinates, density of the densest particle, and number of particles in 
each halo. 

Further analysis is performed by Ganyl, Gadget Analysis Code (Abel et al. 2001) and P-GroupFinder (Springel 2000). 
Ganyl analyzes spherical profiles for GADGET data, taking both the list of centers from HOP analysis and a snapshot of the 
particle data that matches it. It converts units from comoving to proper and sorts according to their gadget IDs. Then, it 
identifies the centers of mass in the snapshot, searches in expanding spherical shells until it finds the virial radius r2oo, and 
thereafter calculates the total halo mass m2oo, velocity dispersion, and angular momentum. 

P-GroupFinder uses a different approach to find halos. The particle distribution is segmented into groups using the 
friends-of-friends (FOF) algorithm. A fraction of the mean interparticle separation is used as the linking length, and halos 
are defined with distance criteria only. Gravitationally self-bound substructures are then extracted from each FOF group by 
using the SUBFIND algorithm. The bound part of FOF halo counts as a subhalo as well, so every FOF halo has its bound 
counterpart. 

The advantage of Ganyl is in its cosmological approach of defining halos using over-densities, but since Ganyl uses density 
centers from the HOP algorithm, it finds over-densities more than once for the same object thereby creating artificial halos. 
Since P-GroupFinder checks for the gravitationally bound structures, it can isolate artificial halos as real world objects. Virial 
radii in kpc are calculated from the total mass in particles and their values are slightly larger than virial radii calculated by 
Ganyl. We use P-GroupFinder to identify black holes as the most-bound particles in their host halos. 

We construct a merger tree from all 33 snapshots. The most massive halo incorporates hundreds of subhalos inside its 
virial radius. We examine all mergers and select a sample of the most massive halos from each group of mergers together with 
a sample of isolated halos that do not merge. For this new list we locate the GADGET IDs for the particles at the centers of 
every halo, and, by tracing IDs in different snapshots, we identify coordinates and velocities to get trajectories. 

Density plots were made with code provided by Naoki Yoshida that applies a spline SPH kernel to derived densities and 
then projects them along a specified axis to produce a two dimensional smoothed image. 

3.3 Black Hole Trajectories 

After identifying a large dark matter halo and its associated substructure, we can track the dynamical behaviour of the 
purported IMBHs in each halo. We follow the evolution of the IMBH population to estimate their final present day distribution. 

We identify the most massive halo in our simulation (SIM1 from now on) at z=l. This halo (the "primary") has a virial 
radius R„ ir =370 kpc and a virial mass M vir = 1.52 x 10 12 M Q . We tracked the primary and its substructure as it assembled 
and tracked the mass and trajectories of the seed black holes. 

If IMBHs are Pop III remnants, then by z > 8, enough time should have passed to allow an IMBH to sink to the center 
of every halo. From the simulation point of view this population can be selected by choosing the most-bound particle in all 
the halos identified by P-GroupFinder. There were 2869 dark matter halos with masses in range 10 7 M Q <5 M <J 10 10 Mq, 
and we label the most bound particle in each halo as an IMBH. By connecting particles' coordinates through 33 snapshots, 
we obtain IMBHs trajectories from z=8.16 to z=1.00. From 2869 IMBHs at z=8.16, 1958 of them can be found inside the 
primary halo at z=1.00. Fig. la shows a density plot of XY-projection of the SIM1 box at z=8.16. For better viewing of the 
Fig. 1, zoom 10 times on online version of the figure. Density peaks (in yellow) are the centers of dark matter halos selected 
by P-GroupFinder (labeled as IMBH). The trajectories of the IMBHs for the most massive 100 halos as they spiral inside the 
primary halo, is overplotted. It is suggested (Miralda-Escude & Gould 2000) that clusters of IMBHs might exist inside the 
inner kpc of every galaxy, spiraling towards the SMBH at the center. Our results are consistent with that point of view. 
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3.4 Black Hole Kick Velocity 

As mentioned in the introduction, coalescence accompanied by gravitational radiation may give an IMBH a significant kick 
velocity. To explore the effect of such kicks on the underlying distribution of IMBH masses, velocities, and occupation fraction 
within forming galaxies, we model these kicks using existing estimates in current literature (Favata et al. 2004). We model the 
distribution of kicks with a truncated Gaussian in the interval {0, 150} kms -1 , with a mean kick of 75 kms -1 . Recall that 
by IMBH, we are referring to the 2869 particles which are located at the centers of sub-halos at z = 8.16; 3/4 of which are 
destined to eventually merge with the primary halo in the absence of kicks. This sets an upper limit to the number of IMBHs 
within this representative large dark matter halo. 

We add this presumed kick velocity, with a random direction, to the velocity of the IMBH at z = 8.16. In this way, we 
create a new snapshot file with changed velocities for the group of black holes only. We use this snapshot file as an initial 
condition file for a new simulation (SIM2a from now on) that starts at z = 8.16 to examine the differences in trajectories of 
the IMBH particle, with and without kick velocities. The same procedure was repeated for a different kick velocity interval 
{125, 275} kms -1 , with a mean kick of 200 kms -1 which was used for simulation - SIM2b. These are the most likely ranges 
suggested in related works (Favata et al. 2004, Merritt et al. 2004). 

Notice the fundamental difference between SIM1 and SIM2. In the case without kicks, black holes embedded in subhalos 
reach the center of the primary halo as dynamical friction siphons energy from the subhalo orbit. The main contribution to 
this process comes from the total mass of their halos which remain mostly bound throughout our simulation. In SIM2, IMBHs 
are generally ejected from their subhalos. Dynamical friction is not efficient on a black hole once it is stripped from a halo, 
due to its much lower mass (Hansen & Milosavljevic 2003, Gerhard 2001, Portegies Zwart 2003). Fig. lb and Fig. lc show the 
density plot of XY-projection of SIM2a box (middle) and SIM2b box (bottom) at redshift z=8.16. Once again the trajectories 
(from z=8.16 to z—1) of a sample of IMBHs are over-plotted. This same population of IMBHs in these three scenarios differ 
only in the IMBHs' velocities and has different trajectories depending on the assigned kick. In the first place, kicks have to 
be large enough to provide IMBHs with a velocity larger than the escape velocity of the host halo. If not, an IMBH will 
not escape the gravitational potential of the host halo leading to a scenario quite similar to SIM1 (compare Fig. la and Fig. 
lb). Larger kicks, as in SIM2b, change the IMBHs' trajectories dramatically (compare Fig. la and Fig. lc). Fig. 2 explains 
the difference. Marked as pluses are the maximum escape velocities calculated from the gravitational potential of every dark 
matter halo selected at z=8.16. Escape velocity decreases with halo radius and reaches its maximum at the center. Hence, 
the maximum escape velocity is the escape velocity of the IMBH since the IMBH is the particle at the halo center. Escape 
velocity increases with the halo mass. Shown as circles are the velocities of IMBHs relative to their host halos at z=8.16 
after assigning them kicks centered at 75 kms -1 . The plot shows that IMBHs from less massive host halos have velocities 
larger than the corresponding escape velocities. With the increase of host halo mass, IMBHs need to acquire larger velocities 
in order to escape from the host halo gravitational potential. For some of them the assigned kicks are not large enough (all 
the IMBHs that lie below the data set represented with pluses). These IMBHs fall back and continue the host halo original 
path toward the primary halo (as when they receive no kicks). This process is demonstrated by the IMBHs in SIM2a (Fig. 2 
circles), hosted by the most massive halos (kick velocities less than escape velocities). Since we plotted trajectories in Fig. 1 
only for IMBHs hosted by the hundred most massive halos, trajectories in Fig. la and Fig. lb arc almost identical. For this 
same set of IMBHs, but now with kick centered at 200 kms -1 in SIM2b (diamonds), almost all have enough velocity to escape 
host halos. As a result, trajectories in Fig. la substantially differ from the trajectories in Fig. lc. Unlike Fig. la and Fig. lb, 
the IMBHs in Fig. lc leave their host halos at z=8.16 and as a result do not form binaries on their path to the primary. 

Even though a substantial kick is assigned to the IMBHs, and dynamical friction does not play an important role for 
them anymore, they still manage to merge with the primary halo. Some fraction of them sink to the center to coalesce with 
the SMBH at the center of the halo. As the gravitational potential of dark matter halos increases at lower redshifts (Fig. 
3), kicked IMBHs are captured and sink to the primary. Fig. 3 shows the maximum escape velocity as a function of dark 
matter halos' mass at different redshifts for all the halos identified in our simulation. The primary halo (thick line) goes 
through a large increase in gravitational potential at lower redshifts. The maximum escape velocity for the primary halo 
increases from ~250kms -1 at z=8.16 to ~700kms -1 at z=l. Hence, the growth of the primary halo and the deepening of its 
gravitational potential well is responsible for the capture of IMBHs originally ejected from their host halos. However, the final 
spatial distribution of the IMBHs is quite distinct (Fig. 4), and their subsequent rate of coalescence with the central SMBH 
consequently may be different by orders of magnitude depending on kics. 

The number of IMBHs inside the primary halo at z=l is 1958 and in SIM2a the number is only smaller by 14 IMBHs. This 
is a coincidence, caused by the fact that we changed the velocities of all the IMBHs identified at z=8.16 and not just of those 
identified inside the primary halo at z=l. This means that some fraction of the IMBHs that did not reach the primary halo 
originally, reached it in this new simulation because of the changed velocities. Clearly the new IMBH arrive "naked", without 
being enveloped in the dark matter sub-halos in which they formed originally. The number of IMBHs inside the primary halo 
in both SIM1 and SIM2a originating from the same sub-halos, is 1851. Thus, 5.4% of the original population found in the 
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primary halo will not be found there due to a change in their velocities at z=8.16. Similarly, a comparable number of the 
IMBHs that originally did not reach the primary halo, will do so in the new simulation. 

The change in the IMBH trajectories in the simulation with larger kicks (SIM2b) is more pronounced (Fig. lc). With 
larger kicks, more IMBHs attain velocities large enough to leave their host halos. Although they need more time to sink into 
the potential well of the primary halo, 1795 of them reach the primary halo at z=l. The number of IMBHs originating from 
the same subhalos as in SIM1 is 1630. Thus, 16.7% of the original population is lost due to the kicks, but for the same reason, 
a new population of IMBHs is kicked into the primary halo, which in total gives 1795 IMBH; implying that the primary halo 
in SIM2b will have only 8.3% fewer IMBHs than the primary halo in SIM1. 



4 POST-MERGER EVOLUTION 

We find that significantly more black holes sink to the center of the primary halo when they are embedded in their dark 
matter sub-halos. However, we find that even large kicks cannot prevent the accumulation of IMBHs at the center of a large 
dark matter halo. Fig. 4 shows the number density of the IMBHs inside the primary halo from SIM1, SIM2a and SIM2b. 
Although the total number of black holes in the primary halo differ from SIM1 by 0.7% for SIM2a and 8.3% for SIM2b, there 
is a decrease in the number of IMBHs in SIM2 for the inner 10% (in radius) of the primary halo. 

We now estimate the dynamical friction time scale for the IMBH and subhalo sinking in the gravitational potential of the 
halo. From Eqn (1), for a singular isothermal sphere with circular velocity v c , the velocity dispersion is a—v c /\/2\ the mass 
inside the virial radius is 1.52 xlO 12 M . For the primary halo these values are r„i r =370kpc comoving, velocity dispersion 
a — 157 kms -1 , and impact parameter b maa; =r„; r . We calculate a radius (i s i„k) which IMBHs have to reach in order to 
merge to the center in less than one Hubble time. We distinguish two cases. First, when the IMBHs are inside sub-halos with 
a minimum mass of 10 7 Mq (SIM1). Second, when an IMBH has been ejected from its parent sub-halo (SIM2). In the first 
case, the IMBH is brought to some radius inside the primary halo by the parent sub-halo and will continue sinking toward 
the center while being embedded within the sub-halo. This will occur for the IMBHs at less then r sin fc=r„; r /30 when the halo 
collapse virialises. We find that 4% of the IMBHs formed are at radii less then r s i„fc. So for this model we predict that, in the 
absence of kicks, 83 IMBHs reach the center and coalesce with the central SMBH (or the seed SMBH formed in the subhalo 
that became the center of the primary halo). More generally, this predicts O(10 2 ) IMBH mergers per Milky Way- like halo 
over one Hubble time; even for the cases shown here where only halos with masses > 2.83 x 10 7 M were allowed to form 
black holes. Since there are ~10 10 galaxies in the observable universe in this mass range, LISA merger rates will be ~10 12 per 
Hubble time or R~100 per year. 

In the second case, the kicked IMBHs have a significantly flatter spatial distribution, partly because they get decoupled 
from their parent subhalos, so there are fewer inside r V i r /30. Fig. 4 shows that there are 2.2% of IMBHs from SIM2a and 0.95% 
of IMBHs from SIM2b inside this radius. Fig. 5 shows the ratio of IMBHs with kicks and IMBHs without kicks. There is a 
large drop in the central population of IMBHs with kicks. This can also be seen in the density plots in Fig. 6. By repeating the 
same calculation for dynamical friction as in the previous case, with the only difference in the value for mass in the equation 
for dynamical friction (8.85xl0 5 Mq instead of 10 7 Mq), we calculate that in order to merge at the center, IMBHs in SIM2 
have to be at radii less then r sin fc=r l ,i r /100. Notice that while we use particles as tracers of the IMBHs, in reality IMBHs 
have smaller masses and need to get to smaller radii for dynamical friction to be efficient. 

Only 0.25% of the SIM2a IMBHs are inside r sin k=r V ir/W0. Thus, under the same assumptions while allowing kicks, only 
about 4 to 6 IMBHs merge with the central SMBH over a Hubble time — a factor of 20 lower merger rate. In SIM2b there 
are no IMBHs inside r„i r /100 except for the one originating from the ancestor of the primary halo at z=8.16. These numbers 
can increase since, in reality, a is lower at small radii resulting in smaller t/ r i c . The merger of 83 IMBHs in SIM1 leads to 
the formation of 7xl0 7 M SMBH inside the dark matter halo with a velocity dispersion a = 157 kms -1 . From cosmological 
Monte Carlo realizations of the merger hierarchy (Volonteri et al. 2003), SMBHs of the same mass and halo velocity dispersion 
can form only when mergers are accompanied with gas accretion as growth mechanisms. However, both the results are below 
empirical values (Ferrarese et al. 2002). We note that the primary halo is not actually spherical and if a halo is triaxial, 
then some fraction of the IMBHs can "walk" into the inner halo (r <JC 10~ 2 r„; r ) region on time scale ~ 10 t or bital due to 
centerphilic box or boxlet orbits (de Zeeuw 1985, Zhao et al 2002, Holley-Bockelmann et al. 2005). In such cases, dynamical 
friction becomes effective in bringing the IMBHs to the halo center. Also, some of the IMBHs are being assigned with kick 
velocities that direct them toward the center. Due to this, some additional IMBHs from SIM1 and SIM2 could reach orbits in 
the center of the primary halo. We will explore the consequences of the non-spherical shape of the primary halo on the long 
term dynamical evolution of the sub-halos, and the implications for black hole mergers and evolution of the inner regions of 
the galaxy in another paper (in preparation). 

Kicks are also responsible for ejecting IMBHs from gas-enriched regions of the halos. Since gas accretion is one of the main 
mechanisms for black hole growth, dumping black holes into regions of lower gas densities would prevent AGN-formation. 
Gas accumulates where the gravitational potential wells are deepest and where dark matter densities are highest. Therefore 
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we may get a crude estimate of the ability of IMBHs to accrete gas by tracing the dark matter density at their positions. Fig. 
7 shows average local density of dark matter traced by our set of IMBHs vs. redshifts. IMBHs get ejected from their host 
halos into the regions with lower density. From z=7.75 to z=2.80, kicked IMBHs are in an environment 1-6 times less dense 
than in the no-kick case. This could suppress formation of AGNs. Note that around z=3.00, the average IMBHs in all three 
cases start tracing similar density distributions. The number density of AGNs at z=8.16 is nAGiv-TMpc -3 assuming that 
every dark matter halo in SIM1 has formed at least one AGN. Ejecting IMBHs in SIM2a results in a lower number density 
of AGNs, nAGJV~10~ 2 Mpc~ 3 (n A Giv~10~ 3 Mpc~ 3 for SIM2b). This implies that for z>3, the contribution of faint AGNs to 
the ionizing background would be decreased if kicks of IMBHs were important. 



5 DISCUSSION 

We have performed and analyzed high-resolution collisionless simulations of the evolution of structure in a ACDM model. 
We have followed the formation and evolution of dark matter halos, and by assuming that each halo hosts an IMBH, we 
have studied the formation of SMBHs. We focussed on two specific cases. First, IMBHs together with their host halos merge 
through dynamical friction. Then, we explored the effect of kicks on the IMBH distribution. We found that kicks lead to 
the ejection of IMBHs from their host halos in many cases. Analytically, it is clear that the dynamical friction will act more 
efficiently on the host halos than on the less massive black holes formed within them. Our illustrative calculations highlight 
some of the expected differences in the final density distribution of the black holes. They can be quite different, even in the 
presence of the modest kick velocities we have imposed. We note that, in general, IMBHs ejected from their parent sub-halos 
at high redshift still get bound to the most massive dark matter halo at later times but they form a "hot population" with a 
flatter density profile than the underlying dark matter. 

In order for dynamical friction to work accurately in numerical simulations, the density of the background (in this case, 
the primary halo density) needs to be resolved well. Higher resolution gives smaller softening lengths for particles and smaller 
softening lengths give more realistic dynamical friction. With the resolution of our current simulations, dynamical friction can 
not be efficiently realized in the model after z~l and r^r„i r /30, because the background density of the parent halos is not 
resolved well enough. However, our simulations are able to account for effects of dynamical friction on the sub-halos at high z 
(Vine & Sigurdsson 1998). Subsequent papers will track the dynamical evolution of the sub-halos and IMBHs at late times, 
using higher resolution simulations and semi-analytic implementations of dynamical friction (Binney & Tremaine 1987). A 
little over 4% of IMBHs merge with the SMBH at the center in less than one Hubble time. If this is the dominant way of 
creating SMBHs, then it is efficient even if the IMBHs are ejected from their parent halos, with merger numbers only reduced 
by a factor of two or so for modest kicks. LISA observations should strongly constrain any kick on IMBHs formed from very 
massive Pop III stars in low mass proto-galactic sub-halos. 

If these IMBHs accrete from the interstellar medium they may also be found by deep X-ray observations (Agol & 
Kamionkowski 2002) in the Milky Way. The observational signature will depend sensitively on the spatial distribution of 
these black holes in our galaxy. Similarly, the next generation of micro-lensing searches may be able to constrain high mass 
black holes from the longest events (Agol et al. 2002). 
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Figure 1. Sample of hundred IMBHs selected from most massive host halos at z=8.16 and their trajectories (white) from z=8.16 to 
z=l overploted on a 2D density projection. Density peaks in yellow correspond to host halos centers at z=8.16 and to the positions 
of the IMBHs. top(Fig la): no kick case (SIM1); middle(fig lb)(SIM2a): case of [0,150] km/s kick centered at 75km/s and bottom(Fig 
lc)(SIM2b): case of [125,275] km/s kick centered at 200km/s. With large kicks IMBHs overcome the host halos' gravitational potential 
resulting in a change of their trajectories and the final distribution of IMBHs (the smearing of the trajectories at the center is quantified 
in the later plots). Observe, e.g., the upper right corner of all three plots where two IMBHs form a binary in SIM1 and SiM2a. In SIM2a, 
IMBHs are trying to escape but the growing dark matter gravitational potential recaptures them. In SIM2b, the assigned kick is larger 
than in SIM2a and also large enough to prevent binary formation. Underneath this binary in SIM1, at the coordinates [5700,7000]kpc, 
a dense halo center hosts an IMBHs. In less then lOOkpc, this halo captures three IMBHs more. This small "cluster" of IMBHs spirals 
in toward the center. The picture clearly resolves IMBHs orbits around center of mass up to the point, [5200,7400]kpc, where the orbital 
separation is smaller then the softening lenght. This point is reached later in SIM2a because of the kicks while in SIM2b binaries never 
form. 
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Figure 2. Black holes' escape velocity as a function of halo mass (z=8.16) as pluses, calculated from gravitational potential of host halo 
set consisting of 2869 members. Black holes' velocity relative to their host halos' velocity reduce to assigned kicks: [0,150] km/s centered 
at 75 km/s repesented as circles and [125,275] km/s centered at 200 km/s represented as diamonds. Lower mass host halos have lower 
maximum escape velocities enabling their IMBHs to escape. IMBHs in more massive host halos demand larger kicks in order to escape. 
As a result, some of them stay captured in their host halo gravitational potential. 




Figure 3. Maximum escape velocity (corresponding to the central gravitational potential) of dark matter halos at different redshifts as a 
function of dark matter halos' masses. The gravitational potential of the primary halo (thick line) increases at lower redshifts, capturing 
even the IMBHs which have been assigned with the highest kick velocities. 
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Figure 4. Number of black holes as a function of primary halo radius. No kick case (thick); [0,150] km/s kick centered at 75 km/s (dash) 
and [125,275] km/s kick centered at 200 km/s (dash-dot). Although the number of IMBHs entering primary halo is comparable in all 
three cases, the difference in their interior distribution is well pronounced. 




Figure 5. Fraction of SIM1 black holes in SIM2a and SIM2b as a function of radius. AN is N - Nkick where N is number of black holes 
in SIM1 (no kick). Thick line for [0,150] km/s kicks centered at 75 km/s and dots for [125,275] km/s kicks centered at 200 km/s. 
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Figure 6. Density of the primary halo at z=l (thick dash) and its most massive progenitor from z=8.16 (dots) as a function of radius. 
Also, density in hosted black holes for no kick (thick); [0,150] km/s kick centered at 75 km/s (dash) and [125,275] km/s kick centered at 
200 km/s (dash-dot). 
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Figure 7. Local density of dark matter traced by black holes as a function of redshift. No kick (thick); [0,150] km/s kick centered at 75 
km/s (dash) and [125,275] km/s kick centered at 200 km/s (dash-dot). Ejection of IMBHs from gas enriched regions of galaxy infuences 
AGNs formation rates, reduces their numbers and their contribution to the ionizing background. 



